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We construct a family of explicit rotational solutions to the nonlinear governing equa- 
tions for water waves, describing edge waves propagating over a plane-sloping beach. 
A detailed analysis of the edge wave dynamics and of the run-up pattern is made 
possible by the use of the Lagrangian approach to the water motion. A graphical 
representation of the edge wave is also presented. 



1 Introduction 

Standing on a gently sloping straight beach, it is a matter of observation that 
various waveforms propagate on the surface of the sea. Among these we find edge 
waves - water waves that progress along the shoreline. These waves, often difficult to 
visualize (this has, probably, prevented the regarding of this waveform as important 
for a long time), are coastal trapped, i.e. their amplitude is maximal at the shoreline 
and decays rapidly offshore. They produce on the beach beautiful run-up patterns 
(highest points reached by a wave on the beach). Although propagation is along the 
straight shoreline and the waveform is sinusoidal in the longshore, these waves are not 
one-dimensional (Lighthill 1978). 

While they were originally considered to be a curiosity (Lamb 1932), edge waves 
are now recognized to play a significant role in nearshore hydrodynamics. For shallow 
beaches empirical evidence shows that incident storm waves loose most of their energy 
trough wave breaking by the time they reach the shore. After breaking offshore, as 
the waves progress to shallower water, their height decreases reaching its least value at 
the shoreline. Since storms often result in pronounced shoreline erosion, the surf zone 
water processes with the onset of a storm are dominated by wave conditions other than 
the incident waves - role attributed by oceanographers to the edge waves ; see Komar 
(1998). There are other instances when edge waves are of significance. For example, 
processed data from the water waves created by an earthquake occuring in April 1992 
in the ocean fioor near the Californian coast show that two distinct wave packets 
(both directly generated in the nearshore by the vertical motion of the ocean bottom) 
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were recorded at a coastal station about 150 km far from the epicentre. At first, less 
than an hour after the occurence of the earthquake, relatively fast-moving swell with 
amplitude around 15 cm striked from offshore. About two hours after the swell has 
subsided, relatively slow-moving edge waves with amplitudes around 50 cm (Boss et 
al. 1995) arrived. Measurements performed on this occasion confirm the rapid decay 
of the amphtude of the edge waves: at an offshore distance of 12 km the amphtude is 
reduced to 10 % of its maximal value (attained at the shoreline). Let us also mention 
that it has been observed (see Evans 1988) that hurricanes travelling approximately 
parallel to a nearby coastline sometimes give rise to edge waves. Interestingly, edge 
waves can in fact be generated directly in a laboratory wave tank (Yeh 1985). 

The edge wave phenomenon has been extensively studied and discussed in the 
mathematical literature within the framework of linear theory. Due to the small 
displacements associated with these waves, the governing equations for water waves 
or the shallow-water equations are linearized (Minzoni & Whitham 1977 showed that 
both approximations are equally consistent) and this simplification permits a thorough 
analysis. We refer to Ehrenmark (1998) for an excellent up-to-date survey. 

Despite the fact that the linearizing approximation lacks rigorous mathematical 
justification, it has been used with considerable success as a large variety of theoret- 
ical studies are confirmed in experimental contexts. The investigation of nonlinear 
edge waves can be seen as a natural extension to the linear theory. Whitham (1976) 
showed the existence of irrotational weakly nonlinear edge waves that propagate par- 
allel to the shore using a formal Fourier series expansion for the full water-wave theory. 
A study of properties of nonlinear progressive edge waves based on the fact that the 
evolution is described by the nonlinear Schrodinger equation was performed by Yeh 
(1985). This paper describes an alternative approach; the main impetus for the results 
reported here comes from the belief that the need for a more rigorous theory remains 
thoroughly justified. A quest for an explicit edge wave solution for the governing 
equations appears to be of interest since the structure of the edge waves obscures 
their visual observation; moreover, a solution in closed mathematical form provides a 
background against which certain features which have been observed (and predicted) 
can be checked. It turns out that the deep water wave solution discovered by Gerstner 
(1809), can be adapted to construct edge waves propagating along a plane-sloping 
beach. This possibility was pointed out by Yih (1966) but the treatement therein, in 
essence followed also by Mollo-Christensen (1982), provides only an implicit form for 
the free water surface. We present a procedure by which exact edge wave solutions to 
the full water-wave equations are obtained. The closed form of the solution in Lagra- 
gian (material) coordinates permits us to provide clear illustrations of the structure of 
these edge waves. From an examination of the solution we also obtain the run-up pat- 
tern, an attractive feature being the occurence of cusps - these shapes are confirmed in 
both field and laboratory evidence (see Komar 1998). Thus, we establish with rigour 
the existence of rotational nonlinear edge waves, unraveling the detail structure of the 
wave pattern. 
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2 The edge wave 



We take a plane beach and adopt a coordinate system as shown below, with the 
shoreline the x-axis, and the still sea in the region 

R — {{x, y,z) : X & W, y < bo, < z < {bo — y) tan a} 

71" 

for some bo < 0; here a e (0, — ) defines the uniform slope. 
Cross-section of the still sea. 



g 




Let u — {ui,U2,U3) be the velocity field and let us recall (e.g. Crapper 1984) the 
governing equations for the propagation of gravity water waves when ignoring viscous 
effects. Homogeneity (constant density p) is a good approximation for water (see the 
numerical data in Lighthill 1978) so that we have the equation of mass conservation 
in the form 

dui du2 dus 



dx dy dz 
The equation of motion is Euler's equation, 



= 0. 



(2) 



Dui 
~Dt 

DU2 

Dt 
Dt 



1 dP 
p dx ' 

1 dP 
p dy 

1 dP 

g cos a. 

p oz 



g sin a. 



where P{t, x, y, z) denotes the pressure, g is the gravitational acceleration constant and 

^/^ n n , Df Of Of df df 

D/Dt IS the material time derivative, —— = — + ui — — h U2 t; — h ^3 — , expressing 

Dt at ox dy dz 

the rate of change of the quantity / associated with the same fiuid particle as it moves 
about. The boundary conditions which select the water-wave problem from all other 
possible solutions of the equations (l)-(2) are (see Johnson 1997): 
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(i) the dynamic boundary condition P = Pq at the free surface, where Pq is the 

constant atmospheric pressure, decouples the motion of the air from that of the water; 

(ii) the kinematic boundary condition at the free surface expresses the fact that 
the same particles always form the free water surface; 

(iii) the kinematic boundary condition at the bottom, requiring the normal velocity 
component at the bed to be zero so that it is impossible for water to penetrate. 

The general description of the propagation of a water wave is encompassed by the 
equations (l)-(2) and the three boundary conditions (i)-(iii), a distinctive feature being 
that the free surface is not known and must be determined as part of the solution. 

We adopt the Lagrangian point of view by following the evolution of individual 
water particles. We suppose that the position of a particle at time t is given by 



where A; > is fixed. It should be pointed out that the quantities a, b, c do not stand 
for the initial coordinates of a particle, they are simply labeling variables serving to 
identify a particle. We may think of them as parameters which fix the position of a 
particular particle before the passage of the wave (in still water) , despite the fact that 
the wave is not developing from the still state - otherwise, the flow would be irrotational 
in view of Helmholtz's theorem (see Johnson 1997) but its vorticity is nonzero (see the 
last section). Let us explain the origin of (3). Gerstner (1809; see Constantin 2001) 
gave the only known nontrivial explicit solution to the full water-wave equations, 
showing that the two-dimensional particle motion (a e R, 6 < 6o < 0, /c > 0) 



represents waves of finite amplitude in water of infinite depth. This suggests that 
it might be possible to construct an edge wave using an approach similar to the 
one for a Gerstner wave field. While the theoretical correctness of this conclusion 
was established by Yih (1966) and MoUo-Christensen (1982), the outcome in both 
trcatcments was an implicit form the water's free surface which makes the obtained 
waveform graphically and computationally inaccessible. The closed form (3) provides 
the full details of the edge wave motion without considerable labour. 

Our aim is to prove that the motion (3) is dynamically possible and that we can 
associate to it an expression for the hydrodynamical pressure P such that the governing 
equations and boundary conditions are all satisfied. The resulting free surface of the 
water will be the edge wave we are looking for. 

The map (3) is a diffeomorphism from the still water region R to the water re- 
gion bounded below by the rigid bed {z — 0} and above by the free water surface 



(3) 
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parametrized by 



_ fcb(l+tana)— fcfeotanc 



(4) { 



" sin (jia + •>/ gk sin a , 
y = 6(1 + tan a) - 60 tan a + ^ gfcb(i+tana)-fcbo tana J^^^ ^ V^A; sin a , 



{bo — 6) (1 + tana) tana + 



tana 
2k 



1 



^2fc(6-6o) (1+tana) 



with a & M., b < bo, and t > 0. Indeed, observing that (a, 6, c) 1— > (a, 6 — c, c) defines a 
diffeomorphism of R^, it is enough to show that the map 




a — — e'^''' sin (^ka + ^/gk~^na 
6' + ^ e^**' cos ^/ca + ^ gk sin a t j 



v 



c + c tan a + 



tana 
~2k~ 



,-2fec(l+cota) 



is a diffeomorphism on R x R_ x R^. To sec this, observe that the third coordinate 
depends only on c, being an increasing function /(c) with /(O) = and hmc_^oo /(c) = 
00. Therefore, we may shce M x R_ x M_|_ by planes parallel to the plane c = and 
in each such plane the particle motion is precisely that of a Gerstner wave field (see 
Constantin 2001). This proves that (3) is a diffeomorphism and it is easy to identify 
the boundary of the image of the region R under it. 

The Lagrangian form of the equation of continuity (that is, the volume-preserving 
property of the flow) is fulfilled since the value of the Jacobian of the map (3) is 
independent of time. This, together with the previously proved fact that (3) defines at 
any fixed time a diffeomorphism, shows that the motion described by (3) is dynamically 
possible. To complete the proof that (3) describes the water motion induced by a 
gravity wave, we have to check Euler's equation (2) and the boundary conditions 
(i)-(iii) for a suitably defined value of the hydrodynamical pressure. 

The acceleration of a particular water particle is 

= (^g sin a e'^^^"''^ sin(A;a + •>/ gk sin a t) , —g sin a e*^^^"''^ cos(A;a + •>/ gk sin a t) , 0^ 
so that the equation of motion (2) is 

— — = —pg sin a e^^^~'^^ sin(A;a + a/ gk sin a t), 

ox 

dP 



—— = pg sin a e^^^ '^^ cos(A;a + a/ gk sin a t) — pg sin a, 

dy 



dP 

dz 



—pg cos a. 
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Passing to Lagrangian coordinates, we obtain the system 
- 

da 

^— — pgsm.a e'^^^^~'^'' — p o sin a, 
do 

BP 

^ — = -pg sin a e^''^'"-"^ - pg cosa + p g cosa {1 + tana) e'''''°e-^'''^^+"''"\ 
with the solution 

P = Po + ^^e^'^'-^^-pg(ccosa+{b-bo) sina) - e-M^^--)e^^^o. 

The obtained hydrodynamical pressure has the same value for any given particle as 
it moves about. At the free surface c = {bo — 6) tan a we have P = Pq so that the 
dynamic boundary condition (i) is satisfied. The kinematic boundary condition at 
the free surface, (ii), holds as at any instant the free surface (4) is the image under 
(3) of the still water surface {c = (6o — b) tana : b < bo}. That there is no velocity 
normal to the sloping shore - this takes care of the boundary condition (ii) - is obvious, 
because at z = we have c = and the motion (3) is planar, without any velocity 
component in the direction of z. The proof of the fact that (3) is an explicit solution 
to the governing equations for water waves on a plane-sloping beach is complete. 



3 Discussion 

We have presented an exact edge wave solution to the full water-wave problem, 
the graphical depiction of which is a fairly easy exercise. Let us now emphasize some 
of its significant properties. We present some simple observations which will provide 
a comprehensive description of this nonlinear wave and the particle motion it induces 
below the water surface. 

The wavelength in the longshore direction A = 27: /k is related to the wave fre- 
quency cu by 

a;^ = gk sin a, 

while the wave period is 
We easily infer that 



T 2x 



\^gk sin a 



A = sm a, 

2n 

so that the length of the edge wave is strongly dependent on its period and to a smaller 
degree on the beach slope. The phase velocity U of the edge wave (4) is given by 
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a fact consistent with the observation that if the bottom is fiat {a — 0), then U — 

and no edge wave exists. The dispersion relation (5) for edge waves is obtained (see 
Johnson 1997) within the confines of the formal linear approximation to the governing 
equations, but in our case the relation is derived rigorously as a byproduct of (3). 
From (3) it is clear that any water particle describes circles as the edge wave passes 

- all these circles lie in planes parallel to the sloped bottom. The radius — e*^*-^"^-* of 

k 

the circle described counterclockwise by a particle is maximal for the particles at the 
shoreline (that is, for 6 = c = 0). 

As pointed out in the previous section, the motion of the water body induced by 
the passage of the edge wave (4) is rotational. The vorticity of the water flow defined 
by (3) is 

^ / du2 dui du2 dui\ 
\ dz ^ dz ^ dx 9y / ' 

by the vanishing of u^. Computing the inverse of the Jacobian matrix of the diffeo- 
morphism (3) as 

/ l+exp[fc(fo— c)] cos k{a+Ut) exp[fc(fe— c)] sin k{a+Ut) 

/ l-exp[2fc(6-c)] l-exp[2fe(6-c)] ^ 

exp[A;(6— c)] sin k(a+Ut) 1— exp[fc(fe— c)] cos k{a+Ut) r\ 

l-exp[2fc(6-c)] l-exp[2fc(fe-c)] ^ 

V i i / 

\ (1+tana) (l+exp[2fe(6o— c— ccota)]) (l+tana) (l+exp[2fe(6o— c— ccota)]) / 

a straightforward calculation yields the expression of the vorticity 

2kU 



at a particle whose parameters are (a, 6, c). Therefore the vorticity is in the opposite 
sense to the revolution of the particles in their circular orbits, decreasing rapidly with 
distance from the shoreline/bed. Observe that, despite the fact that the flow (3) is 
not two-dimensional, the vorticity of each individual water particle is conserved as the 
particle moves about. 

The run-up pattern is obtained by setting 2; = in (4); this forces h — so that 
we have ^ 

— e^^° sin (^a + \J gk sin a , 



(6) 



X = a — — e 
k 



y = bo + — e^^° cos {^a + ^/g]^'^in^ , 



k 

0, 

with a e M. The above formula represents the parametrization of a smooth trochoid 
(if 60 < 0) or of a cycloid with upward cusps (if 60 = 0); it also explains why we 
imposed the condition 60 < as otherwise we would obtain a self-intersecting curve. 
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beach 




Run-up patterns at a fixed instant, viewed in the (x,y)-plane: 
trochoid (above) and cycloid (below). 

beach 




water 



The cusp formation by some edge waves was clearly demonstrated in wave-tank ex- 
periments and found by field measurements on the ocean beach (see Komar 1998). 

Another aspect of interest is the amplitude of the edge wave. To determine the 
elevation with respect to the reference plane 

= ^ + - y) tana : y < bo}, 

we compute the distance of a point {x, y, z) lying on the free surface (4) to this plane, 

/ , \ . sino! r,,, 
d — z cos a-\- [y — Oq) sm a — e °, 

with the understanding that positive/negative values on the right-hand side mean that 
the point hes above/below the plane. Since (with h < bo) 

(7) d = ^ (2e'=^(^+*^"")-'=*°*^""cOs(A;a + ^/^k^^ t) - e2'=Kl+tana)-2fcbotana\ 

we see that the amplitude of the edge wave decays exponentially away from the shore- 
line (as b — > —oo). The same conclusion is reached by a formal linear approximation 
(see Whitham 1979) and explains why edge waves are called "trapped waves". 

As expected - and ensured by (7), the amplitude of the edge wave varies with the 
parameters a e M, b < bo. Prom (7) we also infer that, at a fixed b < bo, the crests 
and troughs correspond to the maximal/minimal values of cos{ka + ^ gk sin a t). At 
a fixed time i > 0, we obtain the crest curves (with m e Z fixed and b < bo playing 
the role of a running parameter) 

2m7r 1 /— — : 

X — — ; — ^J gk sm a t, 

k k 

y = 6(1 + tanQ;) - ftotana + ^ e'=^('+*^"")-'=^°*^°", 
z^{bo- b){l + tan a) + ^ e^'=^°(l - e2'=(^-''o)(i+t-«)), 
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and the trough curves 



X 



(2m + l)7r 1 

T v fl'^ sma t, 



k 



k 



y = 6(1 + tana) - feotana - ^ e^''(^+*^'^")-'=''°*'^°'^, 

A; 



{bo — b){l + tan a) tan a + 



tana 
2k 



g2fcbon _ g2fc(6-bo)(l+tana)> 



Note that for both the crest and trough curves the value of x (at a given time) is 
fixed: standing at that location and looking towards the sea, these curves, orthogonal 
to the shoreline, are fully visible at certain instants. Indeed, taking into account the 
fact that on a crest/trough curve the deviation from the reference plane is 



+ 



sm a 
~2k~ 



^2 g*^^(-'^"'"*^'^'-*)~'^''o tan a 



^2A;6(l+tan a)—2kbQ tan a 



), h<ho, 



in view of (7), the monotonicity of the right-hand side shows that the deviation be- 
comes smaller (in absolute value) with the distance from the shore. This feature can 
be easily recognized in the graphical representation of the edge wave given below. 



The edge wave viewed from offshore. The sinusoidal longshore structure and the 
exponential offshore decay in amplitude are clearly visible. 
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